Spin Fluctuation in Heavy Fermion CeRu 2 Si2 
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Spin fluctuations of the archetypal heavy-fermion compound CeRu2Si2 have been investigated by 
neutron scattering in an entire irreducible Brillouin zone. The dynamical susceptibility is remarkably 
well described by the self-consistent renormalization (SCR) theory of the spin fluctuation in a 
phenomenological way, proving the effectiveness of the theory. The present analysis using the SCR 
phenomenology has allowed us to determine fourteen exchange constants, which show the long-range 
nature of the Ruderman-Kittel-Kasuya-Yosida interaction. 
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Effects of the strong correlation of d- and /-electron 
systems are exhibited in dual aspects of localized and 
itinerant characters El 01- I n heavy-fermion sys- 
tems, observations with energies larger than a small scale, 
Kondo temperature 7k, show local-moment behavior, 
such as the Curie- Weiss susceptibility, and antiferromag- 
netic correlations 0] . While at lower energy scales / and 
conduction electrons form composite quasiparticles with 
a large mass renormalization m*/m oc C/T by a fac- 
tor of up to a few thousands. This large effective mass 
has been ascribed to the local Kondo effect and to near- 
ness to a quantum critical point (QCP) at T = 0, which 
separates the heavy Fermi-liquid (FL) state from an an- 
tiferromagnctic phase of local moments interacting with 
Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange in- 
teractions. Theoretical treatment of the both localized 
and itinerant characters has been a difficult and central 
issue for heavy fermions Q. In fact, a focus of recent ex- 
perimental []| and theoretical studies on systems close 
to QCP, or non-Fermi-liquid behavior, is directed toward 
revealing the nature of the fixed point, i.e., whether it 
is a spin-density-wave type H, M, M or a locally critical 
quantum phase transition [flTloT lllj . 

Attempts to dealing with spin fluctuations in heavy 
fermions including the both localized and itinerant char- 
acters in a wide region close to QCP have been carried 
out by the self-consistent renormalization (SCR) theory 
of the spin fluctuation ||. This extends the SCR the- 
ory, well established for weak ferromagnets and anti- 
ferromagnets of d-electron systems Q, to /-electron 
heavy-fermions using the same form of the dynamical 
susceptibility x(q, E) in the phenomenological way [1@- 
Experimental results of bulk properties, such as temper- 
ature dependence of C/T, have been successfully ana- 
lyzed using the SCR theory for many heavy fermions 
especially for parama gnet ic states H [2] However, an 
experimental test El M of the SCRtheory using the 



archetypal heavy-fermion CcRu 2 Si 2 El El by means of 
neutron scattering and bulk quantities showed that mea- 
sured x{Q}E>) is only semi-quantitatively consistent with 



C/T, posing a question about the phenomenology. Ex- 
tending this work, the present study is aimed at perform- 
ing a more rigorous test of the SCR theory on CeRu2Si2 
by comprehensive measurements of neutron scattering. 
We have shown that the SCR theory describes x(Q> E) 
remarkably well, and that it also provides a rewarding 
method to determine a number of exchange constants. 

CeRu2Si2, which crystallizes in a body-centered tetrag- 
onal structure [the ThCr2Si2 structure, see Fig. da)], is 
an archetypal heavy-fermion compound with enhanced 
C/T ~ 350 mJ/K 2 mol El El- It shows Kondo behav- 
ior with T K ~ 24 K El 

and remains in a paramagnetic 
state down to the lowest temperature investigated. The 
local moments with strongly uniaxial anisotropy develop 
antiferromagnetic correlations with the energy scale of 
k B T K EH3- For T < T K CeRusSiz exhibits FL be- 
havior, for example p — po oc T 2 16;], with renormalized 
heavy quasiparticles, proved by the de Haas-van Alphen 
(dHvA) effect yjl |20j. An intriguing property of this 
compound is the metamagnetic behavior under a mag- 
netic field Hm = 7.7 Tj showing a sharp crossover of 
magnetic states EHHHJUl- 

Neutron-scattering measurements were performed on 
the triple-axis spectrometers HER and GPTAS at 
JAERI. Typical energy resolutions using final energies of 
Ef = 3.1 and 13.7 meV were 0.1 and 1.0 meV (full width 
at half maximum), respectively, at elastic positions. Four 
single crystals with a total weight of 19 g were grown by 
the Czochralski method. They were aligned together and 
mounted in a He flow cryostat. All the data shown are 
converted to the dynamical susceptibility and corrected 
for the magnetic form-factor and the orientation factor. 
It is scaled to absolute units by comparison with the in- 
tensity of a standard vanadium sample. 

The dynamical susceptibility x(q, E) is assumed, in 
the SCR theory 0, to be described by 



Xiq^y 1 = X l(E)- 1 - J(q) , 



(1) 



where the local dynamical susceptibility x~l{E) = 
%L / ( 1 — i-E' /Fl ) , expressing the local quantum fluctuation 
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by the Kondo effect, is modulated by the intersite ex- 
change interactions J r .r', and J(q) = J2r^o Jrfi exp(iq ■ 
r). In the standard treatment y|, Eq. (Qj is expanded 
around an antiferromagnetic wave- vector Q, which is ap- 
propriate and has been used for weak itinerant antifer- 
romagnets of d-electron systems Q, . In stead of this 
expansion, we directly apply Eq. Q for the present anal- 
ysis. The necessity of the nonexpansion form was sug- 
hea vy^fe 

magnetic wave- vectors 013113111. To test the SCR 
theory, hence, we firstly measured Im[x(q, E)} in a wide 
q- and -B-range at a low temperature and fitted the data 
using Eq. Q with the adjustable parameters. Secondly 
we compared the temperature dependence of observed 
Im[x(q, E)\ and C/T with theoretical predictions. 

A number of constant-Q and -E scans covering an ir- 
reducible Brillouin zone were performed at T = 1.5 K. 
The intensity data were fitted using Eq. QJ, i.e., to 



T = 1.5 K HER k f = 1 .22 50'-40'-open-open 



gested for heavy fermions because of much weaker q dc 
pendence 0, E3 ; e.g. CeRu^Si? has three antiferro- 



Im[ X (q,E)] = X (q)T(q)- 



'z: 2 + r(q) 2 ' (2a) 

X(q) = [XE 1 - J(Q)]- 1 , (2b) 
r(g) = XLT L /x(q) , (2c) 

with the adjustable parameters xl, Tl, and J r fi |2*fl| . 
Figure GJa) shows five constant-Q scans selected from 
similar 77 scans. One can see from this figure that the 
fitted curves of these scans well reproduce the experi- 
ments [2(|, considering that x(q) an d T(q) of the single 
Lorentzian, Eq. (|2ap. are determined globally in q space 
by the parameters. 

Figures ^b) and (c) show observed and calculated in- 
tensity maps of constant- E scans with E = 1 meV, view- 
ing the surface of the irreducible Brillouin zone. These 
figures show a good quality of fitting, where we note that 
there are 663 observed data points in the irreducible zone. 
The obtained parameters by the fit of the constant-Q and 
-E scans arc xl = 0.066 ±0.006 (emu/mol Ce), T L ~ 4.4 
meV = 51 K [25J, and fourteen exchange constants listed 
in Table [I] The number of exchange constants, mathe- 
matically equivalent to that of Fourier components, were 
chosen so as to reproduce the three major antiferro- 
magnetic correlations peaked at Aq , fc 2 PJ1 La i and fc 3 
[2.1 12 1| and minor structures at the Z and N points, but 
not to show fine structures due to experimental errors. 
The requirement of more than ten exchange constants 
with both positive (ferromagnetic) and negative (antifer- 
romagnetic) signs reflects the long-range nature of the 
RKKY interaction. To our knowledge, this is the first 
experimental determination of the exchange constants in 
a paramagnetic hcavy-fermion system, which has been 
prohibited for lack of a systematic method. 

Since we have shown that the phenomenological 
Im[x(q, E)] [Eqs. J2J] excellently describes the observa- 
tion at T — 1.5 K, we proceed to the second step at 
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FIG. 1: (a) Constant-Q scans at several wave vectors re- 
duced to the irreducible Brillouin zone, which is illustrated 
on the right side together with the body-centered tetragonal 
structure. Curves are calculated by the fit to the SCR the- 
ory, i.e., Eqs. @. (b) An intensity map of constant-i5 scans 
taken with E — 1 meV at T = 1.5 K is shown on the sur- 
face of the irreducible Brillouin zone. Antiferromagnetic spin 
correlations show three peaks at wave vectors fei = (0.3, 0, 0), 
k 2 = (0.3,0.3,0), and k 3 = (0,0,0.35). (c) Calculated inten- 
sity map using the fit to the SCR theory. 



TABLE I: Exchange constants J r ,o between magnetic mo- 
ments at r = xa + yb + zc and 0. They are defined by 
Eq. Q and determined by the fit of constant-Q and -E scans 
shown in Fig. Q to the SCR theory [cf. Eqs. @]. A positive 
constant J r ,o represents a ferromagnetic coupling with an ex- 
change energy — Jr.o&r&o between Ising variables (ay = ±1) 
(25|l . Errors of exchange constants are ±0.02. 
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FIG. 2: (Color online) Temperature dependence of constant- 
E scans taken with E = 1 meV along the E line Q = (H01) 
(a) and the A line (10L) (b). Curves are calculated by the 
SCR theory using the fit result at T — 1.5 K and the temper- 
ature dependent xl(T). 




1 2 3 
E (meV) 



10 15 

T(K) 



higher temperatures: to compare observations with the 
SCR theory. As discussed in Ref. ||, the most impor- 
tant temperature dependence of x(<7? E) originates from 
that of one parameter Xh(T), and the temperature de- 
pendence of the other parameters can be neglected in a 
low temperature range of k^T <C IY- Under this as- 
sumption the T dependence of Xl(7 1 ) can be evaluated 
by solving 



E 



dE 



exp(/3E) - 1 



lm[ X (q,E)} =n\ (3) 



FIG. 3: (Color online) Temperature dependence of the wave- 
vector dependent susceptibility x(l) an d uniform suscepti- 
bility x c (from Ref. |l(|) (a), and the energy scale T(q) at 
q = fei, kz and the Z point (b). Curves are the calculation 
of the SCR theory using the fit result at T = 1.5 K and 
the temperature dependent \i,(T). Observed x(<?) an d F(q) 
are determined by fitting of constant-Q scans to the single 
Lorentzian form, Eq. I2al . at each q and T. These constant- 
Q scans and fits are shown in (c) for q = k^. In (d) observed 
specific-heat coefficient C/T (from Refs. [l3lll^l ) is compared 
with the calculation by the SCR theory. 



where E c and fx are the cutoff energy and effective mo- 
ment, respectively. This is a sum rule simply implying 
that the total magnetic scattering crosssection integrated 
in q and E spaces is a constant determined by the mag- 
netic moment of the crystal-field ground-doublet. 

Several constant- E and -Q scans were carried out at 
higher temperatures. In Fig.|21we show constant-^ scans 
along the E and A lines (see Fig. . Curves predicted 
by the SCR theory, i.e., Eqs. J5J) with xl{T), show agree- 
ment with the observation within 30 % in a temperature 
range T < 40 K. The constant-Q scans were fitted to 
the single Lorentzian [Eq. ipajl] with x{q) and F(g) de- 
termined at each q and T. Examples of this fitting are 
shown in Fig. Etc) for q — k^. The resulting x(q) and 
T(q) measured at q = k\, fc 3 , and the Z point are plot- 
ted as a function of temperature in Figs. Eta) and (b), 
where the uniform susceptibility is also shown. The pre- 
dicted curves in these figures are in agreement with the 
observations within 30 % up to T = 40 K. From these 
comparisons shown in Figs. [21 and [21 although there are 
systematic discrepancies to a certain extent, we conclude 
that the overall agreement of the temperature depen- 
dence of the observed lm[x(q, E)\ with the SCR predic- 
tions is fairly good in view of the simple assumption. 



The specific heat is dominated by spin fluctuations at 
low temperatures, and has been used as an important 
quantity in applying the SCR theory 0, 0] . Hence it is 
interesting to compare the observed specific heat 0, 
with the SCR theory (27|. The calculated C/T is com- 
pared with the observations in Fig. Etd), which shows 
agreement within 50 %. Considering that the theory em- 
ploys the same method as that for d-electron systems, 
which is not strictly justified for / electrons, the agree- 
ment is fairly good. However, the SCR theory on C/T 
has to be improved for more quantitative purposes, which 
shoul d ag ree with a FL theory of renormalized quasipar- 
ticles [l9| at T = 0. 

We have shown that the SCR theory [j| provides a 
useful description of the spin fluctuations, i.e., x(q?E), 
for the archetypal heavy fermion CeRu2Si2. It seems that 
this phenomenology will be applicable to a certain class 
of heavy- fermions. The success for CeRu2Si2 is partly 
based on the fact that Eq. is a good approximation 
at low temperatures. In fact, almost the same form of 
equation was derived by a l/d-expansion theory [2^ on 
the periodic Anderson model. At high temperatures the 
SCR theory proposed the simple assumption that only 
one parameter Xh(T) depends on temperature, and used 



4 



the exact sum rule of Eq. ©. In the limit T — > oo, this 
recipe gives the correct Curie- Weiss susceptibility x{o) = 
fj? /[k^T — </(<?)]■ This natural extension may account 
for why we obtain the acceptable agreement up to 40 
K (see Figs. |21 and |2l beyond the original expectation, 
T<r L //c B ^51 K |3. 

The three antiferromagnetic wave- vectors fei, fc2, and 
&3 are determined by J(q) = max [cf. Eq. JJJ]. In the 
SCR theory J(q) is ascribed mainly to the RKKY mecha- 
nism and is weakly temperature dependent. However one 
also might attribute the wave-vectors to nestings of the 
quasiparticle bands. According to Ref. [13 > the nesting 
mechanism can be incorporated in Eq. Q by replacing 
J(q) with J s (q) + Jq{q,E), where the first term repre- 
sents exchange interactions involving intermediate states 
of high-energy excitations (including the RKKY interac- 
tion), and the second term involving low-energy excita- 
tions of the quasiparticles. This second term possesses 
the energy scale /ch^k, and consequently would show 
larger T dependence [28j. We speculate that Jq(q,E) 
is smaller than J s {q) for CeRu2Si2 (in zero magnetic 
field), because nesting wave- vectors close to fei, fc2, and 
&3 cannot be easily seen in the Fermi surfaces of the 
band structure 0], and because T dependence of J(q) 
could be neglected in our analysis. On the other hand, 
the quasiparticle mechanism supposedly brings about fer- 
romagnetic exchange interactions under magnetic fields 
close to the metamagnetic crossover 22j. Therefore, to 
resolve these problems on the quasiparticle contribution, 
improved SCR theories based on microscopic models [111 
are awaited. It will be also interesting to apply them to 
those in which the quasiparticles play essential roles in 
spin fluctuations, e.g. UPt3 [3(J and CeNi 2 Ge2 |3lj . 

Singularities of QCPs in itinerant magnets have been 
thought to be described by the SCR theory 0, Q or 
equivalently by the spin-density- wave type fixed point . 
However a recent study of criticality of a heavy-fermion 
system CeCi^-rcAu^ tuned to a QCP 0] caused con- 
troversy on a possibility of a locally-critical fixed point 
ff5L 1 1 1 1] - A key observation, supporting this fixed point, 
is the E/T scaling form Im\x{Q,E)] = T~ a g(E/T) 
at the antiferromagneic wave vector Q [lC|. In con- 
trast, the SCR theory predicts the E/T 3 / 2 scaling form 
lm[x(Q,E)] = T~ 3 / 2 g(E/T 3 / 2 ), which one can see us- 
ing x(Q)" 1 oc T 3 ' 2 at the QCP 2] and Eqs. ©. For 
the present SCR theory of CeRu2Si2, located slightly off 
the QCP, this characteristic T 3 / 2 dependence appears ap- 
proximately as T(Q = k{) ~ A + BT 3 / 2 in a low temper- 
ature range < T < 10 K. Hence we speculate that it 
may be interesting to accurately measure T dependence 
of T(Q) to determine whether it shows the SCR T 3 / 2 - 
dependence or is closer to T 1 -dependence, which would 
suggest the E/T scaling. 

In conclusion, we have demonstrated that the SCR the- 
ory remarkably well describes the spin fluctuations of the 
archetypal heavy-fermion CeRu2Si2- The analysis using 



magnetic excitation data covering the entire irreducible 
Brillouin zone has enabled us to determine the fourteen 
exchange constants. 
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